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Abstract 

We numerically analyze the distribution of scattering resonance widths in one- and 
quasi-one dimensional tight binding models, in the localized regime. We detect 
and discuss an algebraic decay of the distribution, similar, though not identical, to 
recent theoretical predictions. 



1 Introduction. 



The decay in time of the survival probability inside open quantum systems is a 
nontrivial issue, both in Mesoscopic Physics and in Quantum Chaology. Such de- 
cay is determined by the distribution of resonance widths, which has been studied 
extensively ||. 

Still, the effect of localization on the statistics of scattering resonances is not 
completely understood. In the strongly localized regime, some arguments predict 
an asymptotic t -1 law for the probability decay and a T" 1 behaviour for the dis- 
tribution of resonance widths T; more recently, an analytical theory was developed 
||, which slightly corrects the latter into an average decay T -1 ' 25 . 

In this paper we numerically address this question by investigating the T distri- 
bution in a class of quasi-one dimensional models. Hamiltonians in this class are 
given by Band Random Matrices. As such matrices provide models for quantum 
localization not only in disordered solids but in chaotic Hamiltonian systems, too 
[0], our present results are also relevant to the latter class of problems 0. 

We use a computational scheme based on the Effective Hamiltonian (EH) ap- 
proach: however, at variance with the usual way of implementing EHs, which ne- 
glects their energy dependence, we perform an exact (within the limits of numerical 
accuracy) computation of the distribution. Our method is described in sec. 2. In 
this way we indeed detect a large interval of T in which P(T) decays like T~ a with 
a in the range 1 — 2, depending on the model, and on the localization ratio. Our 
results therefore signal algebraic decay, still somewhat different from the predictions 
of ref. |||. Some other differences appear, in the large- r part of the distribution. 

However, comparing numerical results with theoretical predictions is by no 
means an obvious task. First, the theory of ref. is one for a continuous model, 
while our models are discrete. Second, that theory somehow assumes a certain 
ideal coupling to continuum, different from ours. Finally, the large-size asymp- 
totic regime in this problem has some nontrivial features, which impose caution in 
analyzing finite-sample data. 

In sec. 3 we discuss some general features of the problem, based on general facts 
about Anderson localization, and on elementary mathematical estimates which we 
derive for our class of discrete models, and which are well confirmed by our data. 
Our numerical results are described and discussed in the conclusive sec. 4. 
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2 Models, and Effective Hamiltonians. 



We consider one-dimensional lattice Hamiltonians, which describe a wire coupled 
to one or two perfect leads, in the form: 



The first two terms describe the wire and the leads, respectively; the third term 
describes the coupling between them. We label lattice sites by an integer j, in such 
a way that the wire is described by 1 < j < N. The Hamiltonian in the leads 
describes hopping between sites spaced by a fixed integer b; in the basis of vectors 
\j), it has nonzero matrix elements only between sites i, j which belong in the same 
lead, that is, they are both larger than N or smaller than 1. In that case, 



The wire Hamiltonian is a finite, N x N matrix with nonzero elements only between 
sites within the wire. Finally, the operator couples sites i,j spaced by b, one 
lying in the wire and the other in a lead. The corresponding matrix elements again 
have the form (0). 

We have considered two special cases, namely: 

(i) b = 1, a tridiagonal matrix with unit off-diagonal elements and diagonal 
elements given by independent random variables uniformly distributed in the inter- 
val [— W/2, +W/2]. This is a finite sample of a one-dimensional Anderson model 
coupled to leads on both sides. In the discussion below we shall also make ref- 
erence to the one-sided Anderson model, in which the Hamiltonian matrix is a 
semi-infinite rather than a doubly infinite one. This is equivalent to inserting a 
perfectly reflecting boundary at n = N + 1. 

(ii) N > b > 1, the wire Hamiltonian is a Band Random Matrix (BRM), that is, 
a real symmetric matrix of rank N such that Hi j ^ «-> \i — j\ < b. H^j are 
independent Gaussian variables, with variance 4p for off-diagonal elements and (3 2 
for diagonal ones. We have chosen (3 = 1. As analyzed in [f|], this corresponds to 
the "matching wire" regime discussed by Ekonomou and Soukolis 0. 

In both cases (i),(ii) the Schroedinger equation for free propagation in the leads 
has solutions u m = 7^ e%km , with dispersion law E = 2 cos kb. For any energy value 
E in the interval [—2, 2] there are b different allowed momenta, 



H = H il) + H io) + H iio) 



(1) 




'i,j+b + 5., 



%3-b 



(2) 



k 



arccos ~ 2ir 



(3) 



s 
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with s = ... b — 1. There are b incoming and b outgoing waves, hence the wire 
enforces multichannel scattering. It is important to remark that with the choice (§) 
the velocity = ^-\ H=E is the same in all channels^ 

The S-matrix relates amplitudes of incoming and outgoing plane waves, Il,r 
and Ol,r respectively (L and R stand for left and right): 



S 



\Ir) [on) 

In our representation, the S-matrix is a 2b x 2b energy- dependent matrix S s ,i(E), 
where I, s are channel indexes. 

We compute the scattering matrix from the Lippman-Schwinger equation for 
the scattering states 



— Gq ■ Vu* 



u 



(4) 



where u are the free eigenfunctions, w ± the scattering states and V is the "potential 
a matrix of rank N defined via the formula H — Hq 



V . Go is the free Green 



function Gq 



{E ± ie — Ho) ■ It can be computed by a complex integral: 



(Go km =(n\(E- H ±ie)- 1 \m) | e=0 
The scattering matrix is given by: 



1 r2n p i(m-n)k 
' dk 



2ti Jo E±ie — 2 cos kb 



(5) 



<E = 



H=E 



I dkj 
dH 



< Ui\V\uf > 



(6) 



H=E 



where i, j are channel labels, and 



is the density of states in channel i. 



H=E 



Now, G can be written in block form as 





9°i2 


9°n 






git 


9l 


9h 


9° 3 s. 



In this section we use small letters for blocks (i.e, submatrices) and capital letters 
for full operators. Here g^, g® 3 , and gf 3 are semi- infinite matrices, g® 2 is a N x N 

1 In rcf. Q it was noted that the statistics of conductance fluctuations is not significantly 
different with other choices of the dispersion law. 
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matrix, g® 2 and g® 2 are matrices with N columns and infinitely many rows, g 21 and 
g 23 have N rows and infinitely many columns. As the potential matrix is localized 
within the center block (2,2), we write it in the block form 


v 


The poles of the S-matrix are the complex values of energy for which the N x N 
matrix (J — g 22 v) has no inverse; hence, they are given by the roots of 



det(I - glfv) = 



(7) 



In order to solve this equation we introduce an 'effective Hamiltonian' matrix H e ff 
of rank N. This is a well known construction [§§,[§§, but the exact form of the 
effective hamiltonian for our specific models is not immediately derived from the 
general theory. There are in fact certain slight differencies between our effective 
Hamiltonian and the one used in ||, which are probably due to a different choice 
of the Hamiltonian in the leads. Therefore we shall presently give a complete 
derivation for our specific models. 
We start with the identity 

G = G - G VG = (1 - G V)G. (8) 

where G — (E ±ie — i^T) 1 1 e=0 Multiplying on the left by Gq 1 and on the right by 
G" 1 , we have 

G- 1 = Gn 1 - V. 



The block form of (13) is 



9°n 9°n 9% 



92i 9%i 

931 932 



9°23 
9°33 



-9l2 V 
- 922 V ) 
-9l2V 



9u 9i2 

921 922 
931 932 



913 
923 
933 



where 1 and are identity and zero matrices of dimensions corresponding 
position in the infinite matrix. 

For the center block we have g 22 = (1 — 922 v )922, whence, multiplying 
left by g 22 1 and on the right by g 22 , we obtain 



to their 
; on the 



922 



(9°22 



V). 



(9) 
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In order to compute (#22) 1 we start from the definition, which in block forms reads 

U12 



9°n 9° 12 9°i3 

921 9%2 923 

9si 932 933 



E-h° n 

U21 





E-h° 22 
U32 





U23 
E-h% 





"1 





0" 







1 













1 



where = —h^ are the coupling matrices. 
The equation for the center block is 

021«12 + 92 2 ( E ~ h 22) + #23^32 = 1- 



(10) 



From (^]) we know that the free Green function is a Toeplitz matrix, whose 
center block g 22 can be written: 



x 

XN-l 



Xl 



XN-l 
Xl 

x 



where x s = if s is not a multiple of b, while, if s = lb, then xib = Xoe dk , with 



k = arccos (-f ). Keeping this in mind we get 



(#21 M 12 +023^32 ) 



Xb 



Xl 







x N+ i 



XN+b 



Xb 



-922- K b . 



_x N+b ... x N+ i ... xi 

where K b = diag(e lkb , . . . , e lkb , 0, . . . , 0, e lkb , . . . , e lkb ) is the "self-energy" . 
Then from ([H]) we obtain: 



(9: 



\-l 
22 J 



(e — h 



o 

22 



K h ). 



Finally, (Eh yields: 



(9 



22 - 



(9 



N-l 

22 J 



v = (E-h% i -K h -v) = {E- h 22 - K b (E)) = (E — H eff {E)). 



where H e ff(E) = h 22 + Kb(E) is the effective Hamiltonian matrix of rank N. The 
role of the effective Hamiltonian emerges on noting that: 



I ~ 922 V = 9l2{{9. 



22) 



v) = g° 22 (E - H eff (E)) 
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which shows that solving eqn.(^) is the same as solving the equation: 

det(H eff (E) - E) = (11) 

For further use we rewrite H e ff in operator notation. In place of h 2 2 we rewrite 
H^: the Hamiltonian operator of the sample, with Dirichlet conditions at n = 
and n = N + 1. Then 

( b N ) 

H eff {E)=H^ + e^ b \Y J + E (12) 

[n=l n =N-b+l ) 

The leads only affect diagonal elements, at sites near the contact points. For the 
Anderson model, where hopping only occurs between neighbouring sites, only the 
first and the last diagonal elements of the Hamiltonian are affected. Some straight- 
forward modifications to the above construction are necessary in the one-sided An- 
derson case; we omit details here. 

Solving the nonlinear eqn.([ll]) is a difficult task. When using the effective hamil- 
tonian formalism, one typically neglects the dependence of H e ff on energy, so the 
problem is reduced to finding eigenvalues of the effective hamiltonian at a chosen 
fixed value E of the energy. Such eigenvalues depend on E Q as a parameter; for 
convenience of language, we will term them parametric resonances in the following, 



reserving the name exact resonances to solutions of eqn.(|TT|). 



3 Theoretical Premises. 

A few remarks are in order, about the mathematical problem set by the above 
formalism. These are most simply formulated for the one-sided Anderson case, so 
we restrict to that case; nevertheless, similar arguments can be developed for the 
two-sided Anderson and for the band matrix models. 

If E is a root of eqn. ([TT|) , then there is a vector of unit norm, satisfying the 
equation: 

H®\iJ)) + e ik W\l)(l\,J>) = E\1>) (13) 

where is now the Anderson Hamiltonian with Dirichlet boundary conditions 
at n — 0, n = N + 1. Multiplying eqn. ([13|) on the left by ?/>, and taking imaginary 
parts, we get: 

lm(e lk(E) )\(l\^)\ 2 = T (14) 
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where we have set E = x + iT with real x, T. 

From the dispersion law the following analytical expression of e 1 ^^ follows: 



e , ME) = E + VWrr-4 (i5) 

In the physical sheet, the square root in eqn.(|l^) has to be chosen such that its 
imaginary part is opposite in sign to V — Tm{E)\ then an easy computation shows 
that Im(e ,fc ' E ') is opposite in sign to T, so eqn. flliD cannot have solutions with 
r ^ 0. In order to solve it, we have to continue H e ff into the nonphysical sheet, 
across the branch cut [—2, +2]. 

We now address the problem of finding estimates for the largest T. We consider 
parametric resonances first: if E is chosen in (—2,2), then eqn.([TJp immediately 
sets a sharp, realization-independent bound on parametric V: 



A-El 



\Y\ < | sm(k(E ))\ = v 2 < 1 (16) 

For exact resonances we can only establish a milder, realization-dependent bound. 
Multiplying eqn. (|i~3D on the left by (1|, using eqn.([T3|), and the explicit form of 

we get: 

m) = (l\iP)(E-e^-V(l)) (17) 



where ^(1) is the random on-site potential at site 1. Since |(2|?/>)| 2 + |(l|-?/>)| 2 < 1, 
using eqn.(^) we get the inequality: 

' < . ._ L (18) 



lm{e ik W) ~ 1 + \E - e ifc ( £ ) - V(l)\ 2 

which, given a realization of the random potential, sets an upper bound to V. At 
the center of the spectrum and small V(l) this bound has approximately the form 
|T| — In turn, this implies that those realizations of the random potential 

which yield T's larger than a given 7 have a probability not larger than ~ I/7. 
This bound on the large-r behaviour of the -P(T) distribution of exact resonances 
is much milder than the bound (|16|) for parametric resonances. We anticipate that 
this difference is manifest in numerical data (see fig. 4, 5 and 6) 

In the limit N = 00, both H® and H eff (E ) become operators in 00— di- 
mensional Hilbert space. The latter operator is obtained by adding a rank-one 
perturbation to the former, which in turn has a dense point spectrum in the inter- 
val I = [-W/2 — 2, +W/2 + 2] (with probability 1). From general operator theory 
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it follows, that in the limit N = oo the distribution of parametric T's collapses into 
a Dirac delta at zero - unlike other scattering statistics, (e.g., the phase shift dis- 
tribution), which have a smooth limit distribution. This physically intuitive result 
should be valid for the distribution of exact resonances, too. 

Thus, on increasing N, we should expect the leftmost part of the finite- N dis- 
tribution -P(r) to rise, and the rest to gradually subside. To further illustrate this 
point, we reformulate eqn.([T3|) by projecting it onto the basis of eigenvectors u n of 
Denoting E n the corresponding eigenvalues and ip n the amplitudes of ip, we 

get: 

e ik ^(u n \l)(l\ip) 



Vv. 



E,„ — E 



whence: 



-W>e 



ik(E) 



\{Un\l)? 
E n — E 



As (ll^) = is excluded, the resonant values of E must solve the key equation: 

E = - e ~ lk{E) ( 19 ) 

where p n = |(l|w n )| 2 . 

At large N, the eigenfunctions exponentially localized, with localization 

lenghts £,{E n ). For resonances E = x + iT with |T| 1/N, a single-pole approx- 
imation to the lhs of eqn.(|l^) should be valid, because eigenfunctions with E^s 
much closer than the average level spacing typically have an exponentially small 
overlap; so the sum in is dominated by a single term|Tl]. Hence, the narrowest 



resonances can be assumed to solve 

Pn e tk ^ &{E- E n ) 

with p n small, so 

E^En+p n e ik{En \ and F w p n sm{k{E n )). (20) 

If we further restrict near the center of the spectrum, then the smallest T come from 
states u n localized around sites n Q (n) lying in the rightmost part of the sample. For 
these, log(p n ) ~ —2no(n)/l; + i] n , with £ the localization length at the band center 
and rj n a fluctuating quantity of order y/n . 

Thus the distribution of very small T's is ruled by the distribution of p n 's. 
If in addition F exp(— 2iV/£), the latter distribution is mainly determined by 
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fluctuations of r] n ; assuming a gaussian distribution for the latter, one gets that in 
this region P(T) has the lognormal distribution already well known in this context. 
Nevertheless, this part becomes negligible at large N, because it comes of a fraction 
~ £/iV of the full set of all resonances. 

Away from this extreme region, the T statistics is more and more affected by 
the change in n . Neglecting r\ n completely, one deduces a dependence ~ 1/T, by a 
simple argument already reproduced in ref. ||. The presence of r\ n just smoothens 
the cusp of 1/T at T pa exp(— 2iV/£), but the 1/T law again re-emerges at larger T. 

So finite-N, normalized distributions P(T) have a peak at ~ exp(— 2N/£), of 
height ~ exp(2iV/£). It is this very peak which eventually builds the limit 5— 
distribution. In the range exp(— 2iV/£) < T <C 1/N, the above rough argument 
suggests a law ~ 1/T; but it must be mentioned that according to Titov and Fyo- 
dorov (H) this inverse law should be restricted to a smaller range. Some numerical 
data about this issue will be given in the next section. 

The large T region is essentially determined by the coupling to continuum, so 
it should be model-dependent. Nevertheless, it is reasonable to assume that the 
number of resonances involved is constant, of order £; if so, again this tail should 
subside at large N, at the rate ~ £/N. 



Finally we use eqn. ([T9|) to investigate the reliability of parametric resonances as 
approximations of exact ones. The equation for parametric resonances is obtained 
from eqn. (p~9f) by replacing in the rhs E by a fixed E chosen in (—2, +2). An 
obvious requirement for parametric resonances to approximate exact ones is that 
their dependence on E be mild. Now, 



dE 




dE 





E 



IWi)p 



(E n - E)< 



dk(E 



dE 



> 



4 



F 2 



If a given parametric resonance is to vary little on changing Eq on the order of its 
width T, it is therefore necessary that: 



r 2 



A-El 



which shows that the parametric approximation becomes unreliable close to the 
edges ±2, and in any case at Y ~ 1. 
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4 Numerical Method, and Results. 



As we are only interested in the statistical distribution of the imaginary parts of 
solutions of (|TTD , we don't need to compute them exactly. Instead, we use the fact 
that the number of zeros of an analytic function f(z) inside a closed path is equal 
to the variation of the phase of the function itself along the path divided by 2n. We 
have therefore considered rectangular regions {E : \Re(E)\ < Eq, —7 < Im(E) < 0} 
in the lower part of the 2nd Riemann sheet. By numerically computing the phase 
of det(H e ff — E) along the boundaries of such regions we obtained iV(7,_E ), the 
number of resonances having real parts in (— Eq, +Eq), and widths not larger than 
7. Typically E = 0.5 in our computations. 

Repeating the procedure for different realizations of our random Hamiltonians, 
we obtained the histograms of resonance widths shown in Figs. 1-6 (there P(T) 
is the probability density for T values). The numerical procedure is quite heavy, 
so we were able to process at most 400 realizations of the BRM model. With the 
Anderson model, computation is faster, so we could process up to 1000 realizations. 

Most of the distributions of resonance widths P{T) computed by the above 
discussed method decay, at very large T, faster than power-like, also because of the 
difficulty of numerically building good statistics in this poorly populated region. 

However, in an intermediate region of values of I, the observed decay is algebraic, 
proportional to T~ a . The width of this region depends on the localization ratio r, 
which is proportional to b 2 /N and to 1/(NW 2 ) for the BRM model and for the 
Anderson model respectively. The region of algebraic decay is very broad in strong 
localized systems, r <C 1; it shrinks as r is increased, and eventually disappears in 
the metallic region r ^> 1. 

The behaviour of the exponent a is somewhat different in the BRM and in 
the Anderson model. In the former case, a increases as r increases, going from 
values about 1.5 — 1.6 to values about 2 — 2.1 (see fig.l). Data obtained at fixed 



r and different values of N, b ( with y > 5, though) show that a only depends 
on the localization ratio in the explored parameter range: see fig. 2, which shows 
distributions with approximately the same localization ratios and different N, b. 
We have numerically computed distributions of widths for the one and two- 



leads Anderson model, too, by solving eqn. (|il|) as explained above. In this case, 
a remains more or less constant around 1.7 — 1.8: it doesn't seem to depend on 
the localization ratio (fig. 3) in the explored parameter range. The only effect of 



10 



decreasing W is the predicted shift of the peak of the distribution towards larger 
values T. 

Most of our data do not yet pertain to the true asymptotic large- N regime. In 
fact, the peak in the left-hand part of our data is still relatively broad in comparison 
to the right-hand part, so the asymptotically interesting region e~ 2N ^ C T C 1 
is still somewhat restricted. Nevertheless the region of algebraic decay is already 
clean, and relatively stable against variation of the localization ratio. We therefore 
surmise, that in more localized regimes the local exponent a would be found to 
smoothly decrease with T, tending to 1 as the peak of the distribution is approached 
from the right. As we shall presently discuss, we get indications in this sense from 
distributions of parametric T's, whose computation could be pushed to significantly 
more localized situations than accessible to exact computations, based on solving 
eqn.(0). 

Parametric distributions were obtained by diagonalizing H e ff(0). Both for the 
BRM and the Anderson model, they exhibit a cut-off at F ~ 1 which is absent in 
the real distributions. The latter in fact decays to zero much more gently. In other 
words: in all the models we have studied, neglecting the energy dependence of ef- 
fective Hamiltonians yields acceptable results only at resonance widths appreciably 
smaller than 1 (figs. 4, 5 and 6). The reason is that, at such large widths, the 
effective hamiltonian is significantly changing already over the width of single res- 
onances. These numerical findings are fully consistent with the analysis in section 
3. 

Exact and parametric distributions fairly well agree in the central part, on the 
right of the peak, but they are again somewhat different on the left. The reason 
is probably that our exact distributions collect resonances from a relatively broad 
interval of real energies, so they include resonances closer to the edges, where smaller 
T's are a priori expected (eqn.(pO|)). 

The sharp cut-off of numerical distributions of parametric resonances was also 
observed in ref. ||, and found to be consistent with the therein developed theory 
for exact resonances in a continuous, white-noise model. That coincidence between 
exact and parametric distributions at relatively large T raises an interesting theo- 
retical problem about the role of the coupling. |L2] 

In fig. 4 we compare parametric and exact distributions both for the one- and 
the two-sided Anderson model. Note that one-lead data refer to a significantly 
more localized regime than two- lead ones, and yield smaller a's. The estimated a 
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for one-sided Anderson is in this case wll In fig. 5 we have again computed the 
parametric distribution for one-sided Anderson, this time in an even more localized 
situation. Moving from left to right, the lhs part of the distribution now exhibits 
a ~ 1 over three decades, after which the distribution gradually drops to zero in 
roughly two decades. 

We finally note that the rightmost part of the exact distribution in fig. 5 is well 
fitted by a T~ 2 law. As the latter coincides with the upper bound established in sec. 3 
for exact resonances, we have an indication that that bound is probably optimal. 

In summary: in this paper we have analyzed the distributions of imaginary parts 
of resonances in different tight-binding models. The most interesting result is the 
presence of a region with power like decay, both for Anderson and BRM model. 
It was obtained by numerically implementing the effective Hamiltonian method, in 
a way which doesn't a priori neglect energy dependence. Our analysis indicates 
that most of the distribution of widths tends to concentrate within a single peak 
around Y ~ exp(— 2iV/£) as more and more localized regimes are approached. While 
subsiding (with increasing N), the distribution of widths on the right of the peak 
displays an algebraic dependence on T. The related exponent ranges from 1 (at 
small T), to 2 (at large T). In the latter region, however, the specific form of the 
coupling to continuum plays a role. 

Useful discussions with G.Maspero are acknowledged. We are grateful to Y.Fyodorov 
and M.Titov for communicating us unpublished details of their work. 
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Figure [I]: Distribution of exact resonance widths for BRM models with different 
localization ratios. Diamonds: N = 100, bandwidth 5, the slope of the dashed 
line is —1.6. Squares: N = 150, bandwidth 26, the slope of the solid line is —2.1. 
Crosses: N = 200, bandwidth 20, slope -1.85. 

Figure [2|: Distribution of exact resonance widths for BRM models with similar 
localization ratios. Diamonds: N = 150, bandwidth 10. Squares: N = 200, 
bandwidth 12. Circles: N = 200, bandwidth 7. Triangles: N = 150, bandwidth 6. 
Their localization ratios S are 0.67, 0.72, 0.245, 0.24, respectively. 
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Figure ^ Distribution of exact resonance widths for the two-leads Anderson model, 
with different localization ratios. Diamonds: N = 150, W = 1.1. Circles: N = 100, 
W = 1.5. Crosses: N = 100, W = 0.3. The slope of the line is -1.75. 
Figure f|: Distribution of parametric and exact resonance widths for the Anderson 
model. Diamonds: exact resonances at iV = 100, W = 2.5, one lead. Squares: 
exact resonances at N = 100, W = 1.5, two leads. The dashed and the solid lines 
represent parametric resonances obtained by diagonalizing H e ff(0), for 1- and 2- 
leads model, respectively. 

Figure |5] Distribution of parametric and exact resonance widths for the one-lead 
Anderson model. Squares: parametric resonances at N = 300, W = 2. The slope 
of the main straight line is —1, that of the shorter one on the right is —2. Diamonds 
are reproduced from Fig.^, which refers to a less localized situation. 
Figure |6|: Distribution of parametric and exact resonance widths for the BRM 
model. Diamonds: iV = 100, bandwidth 15. Circles: iV = 100, bandwidth 5. 
Triangles, N = 150, bandwith 10. The curves represent parametric resonances 
obtained by diagonalizing H e ff(0). 
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Figure 1: P(T) for BRM model 
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Figure 2: P(r) for BRM model 
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